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ABSTRACT 

The gust response of a 2-D cascade is studied by solving the full nonlinear Euler equations employing 
higher order accurate spatial differencing and time stepping techniques. The solutions exhibit the 
exponential decay of the two circumferential mode orders of the cutoff blade passing frequency (BPF) 
tone and propagation of one circumferential mode order at 2BPF, as would be expected for the flow 
configuration considered. Two frequency excitations indicate that the interaction between the frequencies 
and the self interaction contribute to the amplitude of the propagating mode. 


1. INTRODUCTION 

One of the main contributors to fan tone noise is 
the so-called rotor-stator interaction noise. It is the 
result of the interaction of the periodic 
disturbances of the rotor blade (mean) wakes with 
the stator vanes. One method of computing rotor- 
stator interaction was described in detail in [1], in 
which a two-dimensional strip description of the 
unsteady aerodynamic interaction between the 
rotor wakes and the stator vanes is combined with 
a 3-D acoustic response of an annular cascade to 
an incident gust. To improve the computation of 
unsteady aerodynamics, a time-linearized analysis 
method is used in [2], where a linearized Euler 
analysis is shown to produce reasonably well the 
first two harmonics of the blade passing frequency 
(BPF). On the other hand, nonlinear time-domain 
analyses would enable understanding of linear- 
nonlinear regimes, self interactions and 


interactions of two-frequencies. At NASA Glenn 
Research Center, an effort is underway to develop 
a nonlinear time marching algorithm to examine 
rotor wake - stator interaction noise. In the present 
paper a non-linear time marching approach is used 
to compute the gust response of a 2-D cascade. 

The gust cascade interaction problem can be 
analyzed employing either a time-linearized 

approach or a nonlinear time marching approach. 
A number of time-linearized analyses have been 
carried out starting with the Sears problem of an 
airfoil encountering a gust. Some of these studies 
are summarized in [3]. In a time-linearized 
analysis, the gust is assumed to be a small 
harmonic perturbation of the uniform steady mean 
flow, yielding a system of equations in the 
frequency domain. However, in the recent 

linearized Euler analyses (e. g. [2]) fully 

nonuniform mean flows are considered. In the 
time-linearized, frequency domain approach, 
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typically the gust response to a single frequency 
excitation is studied. Nonlinear effects cannot be 
assessed when the amplitude of excitation is no 
longer in the linear range. 

In a time domain analysis, the full nonlinear Euler 
equations are solved using a higher order accurate 
time marching scheme. Of course, the time 
marching technique requires significantly larger 
computational time compared to the linearized 
Euler analysis. But, the advantage of a time 
domain analysis is that from a single solution for 
multi-frequency (harmonics) excitation, responses 
at all harmonics of interest can be extracted. Also, 
nonlinear interactions [4], both self interaction 
generating higher frequency responses and 
interaction between frequencies of excitation 
resulting in sum and difference tones for multi- 
frequency excitation may be studied. Thus, a true 
nonlinear time domain analysis has the distinct 
advantage of being able to mimic the real flow 
more closely. However, to be able to resolve the 
nonlinear effects, higher order accurate spatial 
derivative schemes, time marching schemes, and 
boundary condition treatments are essential. 

The computational time requirements and 
resolution constraints have so far limited the time 
domain analysis to single airfoil gust response 
analyses [5-7] and flat plate cascade [8-9] studies. 
In the present paper, a cascade representing a 
typical modern high-speed configuration and a gust 
typical of the measured wake flow behind a 
modern high-speed fan rotor are considered. 
Following the description of the cascade problem, 
the solution procedure is described briefly. The 
results for single and two-frequency excitations are 
presented and discussed. Nonlinear effects 
observed in single and two-frequency excitations 
are also examined. 


2. THE PROBLEM 

The 2-D cascade geometry shown in Figure 1 is 
the unrolled section at a radial location of a 
modern high speed turbofan stator vane geometry. 
The cascade has a gap-to-chord ratio of d/c = 2/3 
with the inflow and outflow planes located at 
x = 3/2 c, where c is the vane chord, and d is the 


vane gap. The mean flow conditions at the inflow 
and outflow planes are given as: 

Inflow: Pi =1, Ti = 1, and a, = 36° 

Outflow: p 0 / Pi= 0.92 

where Pi and Ti are the normalized mean stagnation 
pressure and stagnation temperature respectively. 
OCi is the mean flow angle and p 0 is the mean static 
pressure at the outflow plane. The inflow periodic 
wake disturbance is described at the inflow plane 
as follows: 

u g (y, t) = {flj cos (k y y-cot) + 
a 2 cos(2 (£ y y - cot) + (p 2 ))}e p 

P g (y,t) = 0, p g (y,t) = 0 

e p = cos (p)e x - sin (J3)e y , j3 = 50° 

co=37i/4, k y = 1 ln/9 , ai = 0.007, 

a 2 = 0.003, cp 2 = -771/5. 

where co is the fundamental reduced frequency, k y 
is the transverse wavenumber, and a/s are the gust 
harmonic amplitudes. The frequency is normalized 
by the chord divided by the ambient speed of 
sound, wavenumber is normalized by the vane 
chord, and gust amplitudes are normalized by the 
ambient speed of sound. The harmonic amplitudes 
of excitation ai (at BPF), and a 2 (at 2BPF) have 
been varied to study the nonlinearities associated 
with the high amplitudes and two-frequency 
excitations. 


3. GOVERNING EQUATIONS 


The full nonlinear Euler equations governing the 
2-D cascade flow are solved in the present 
investigation. In Cartesian coordinates these 
equations are written as: 


dQ dE dF n 
— + — + — = 0 
dt dx dy 


( 1 ) 


where 
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Q = 


P 

pu 

pv 


E 


( 2 ) 


pu 

pu 2 + p 
puv 

u{E + p ) 

pv 

puv 

pv 2 + p 
v(E + p) 


( 3 ) 


( 4 ) 


code solves the flow equations written in chain- 
rule curvilinear form, using a prefactored sixth- 
order compact scheme for spatial differencing, 
with 11 -point explicit stencils at the block 
boundaries. An explicit 10 th order constant 
coefficient artificial dissipation is used for 
damping. The time marching uses a 2N Storage 
fourth-order nonlinear extension of Hu’s 5-6 Low 
Dissipation and Dispersion Runge-Kutta (LDDRK) 
scheme [13,14]. The inflow and outflow boundary 
conditions are implemented using Giles non- 
reflecting boundary conditions [15]. On the 
airfoils, the time derivative of the velocity normal 
to the wall is set to zero. 

In the present computations, grid singularities 
have been removed. However, it was shown in 
[16] that grid singularities have negligible effect 
on the converged solution. The technique of 
implementing mean flow boundary conditions are 
discussed in [17] and the effects of higher order 
wall boundary conditions on the solution are 
explored in [18]. 


and 

P = (y-l){E-^p(u 2 +v 2 )} (5) 


These equations are cast in generalized curvilinear 
coordinates employing the chain rule formulation 
as: 

dr dxd£ dxdr/ dyd£ dydr/ 


4 . SOLUTION PROCEDURE 

The full nonlinear time dependent Euler equations 
governing the 2-D cascade flow are solved 
employing a new parallel computational 
aeroacoustic (CAA) code developed by Hixon 
et al. described in part in [10]. Details of the 
computational scheme may be found in [11-12]. 
Since the flow has a wide range of length and time 
scales, a higher order numerical scheme is needed 
to reduce dissipation and dispersion errors. The 


5. NUMERICAL SOLUTION 

The fan stage considered in this study has a 22 
bladed rotor and a 54 vane stator. The propagating 
acoustic modes are determined from the Tyler- 
Sofrin [19] criterion, m = n N b - k N v , where m is 
the circumferential mode order, n the harmonic 
number, N b number of rotor blades, N v number of 
stator vanes and k a positive integer. For the fan 
stage considered here, the blade passing frequency 
(BPF) tone is cutoff. That is, the BPF interaction 
mode m= -32, and the rotor locked mode m=22 
are expected to decay exponentially. The 
propagating mode is m = -10 at 2BPF. Because of 
periodicity, it is enough to simulate one half of the 
fan. In the rolled out 2-D cascade, only 27 stator 
passages need to be considered and accordingly the 
propagating mode is also stated as m = -5. The 
decaying modes at BPF are m= -16 and m=l 1 for 
this cascade. 

The grid used in the computations has nine grid 
blocks per passage (Figure 2), with 12704 grid 
points, and 27 passages are solved as mentioned 
above. An initial grid tested had 6 grid blocks 
(with 7004 points per passage) which provided 
reasonable flow resolution and acoustic response 
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amplitudes. However, the grid shown in Figure 2 
was chosen to significantly increase the flow 
resolution in the leading edge, trailing edge, and 
wake regions. The grid density varies depending 
on the region of the flow. In this grid, the 
minimum number of grid points in the axial 
direction is 124 and in the lateral direction is 46. 
The minimum number of points per wavelength at 
2BPF is 41 and thus the grid density is sufficient to 
produce the resolutions required for the excitation 
frequencies considered. For each set of gust 
amplitudes, the solution is run till periodicity in 
pressure is achieved on the inflow, outflow 
boundaries and on the airfoil surface. In terms of 
the period of the BPF harmonic, it was found that 
the solution is to be run about 200 periods. Once 
the periodic solution is established, it is post- 
processed to extract the modal information. A 
joint spatial - temporal Fourier transform of the 
acoustic pressure p(x,y,t), provides the amplitude 
of the acoustic waves as a function of rotor 
harmonic n, spatial mode order m, and axial 
location x [20]. To perform a Fast Fourier 
Transform (FFT) the number of samples employed 
is usually a power of two. Since FFT is done over 
the entire computational domain, the number of 
samples cannot be too large. But, it was found that 
a minimum of 32 samples (points in the time 
series) is needed to avoid aliasing problems. 


6. RESULTS AND DISCUSSION 

The time domain analysis is employed to study the 
acoustic response to single frequency and two 
frequency excitations. Single and two-frequency 
excitation responses are compared and linear and 
nonlinear characteristics are examined. 

2BPF excitation : 

The unsteady flow features that are observed in the 
present gust-cascade interaction problem for a 
2 BPF, single frequency excitation is first 
presented in Figure 3. As described above the 
solution is carried out for the 27 passages. 
However, to illustrate the unsteady flow 
characteristics only three full passages are shown. 
The imaginary part of unsteady u-velocity contours 
(Figure 3a) indicate the nature of the gust 


convected through the passage. The wake behind 
the trailing edge is thicker than that would be 
observed in a linearized solution. This is also 
illustrated in the instantaneous u-velocity (Figure 
3b) and vorticity (Figure 3c) contours. The 
oscillations in the vorticity contours from the 
airfoil surface and the wake are typical of high 
accuracy fine grid solutions. Computed results 
indicate that the thickness of the vortex sheet 
behind the trailing edge and the magnitude of the 
vorticty in the wake are functions of the amplitude 
of excitation. 

Figure 3d shows the contours of the imaginary part 
of unsteady pressure. The circumferential mode 
order (m) pattern is observed. The mode order 
m= -5 is propagating both in the inflow and 
outflow regions. The clean and nearly straight 
contours near the inflow (x= -1.5) and near the 
outflow (x=1.5) planes suggest that the reflection 
from the boundaries is minimal. The high 
accuracy spatial differencing scheme and the time 
stepping scheme contribute to the quality of the 
solution. 

The circumferential mode pattern in the entire 
computational domain is shown in Figure 4a. The 
propagating circumferential mode order, m= -5 is 
clearly seen both in the inflow and outflow 
regions. As indicated above, at the inflow and 
outflow planes the reflection appears to be small 
and seems not to affect the solution. The 
propagating mode amplitude as a function of axial 
distance x is shown for the inflow region in Figure 
4b and for the outflow region in Figure 4c. 

The convergence of the pressure waveforms are 
examined for the establishment of the periodic 
solution. The typical pressure waveforms at the 
airfoil pressure and suction surfaces, inflow plane 
and outflow plane are shown in Figure 5. The 
figures show the waveforms at two successive 
periods, 203T and 204T. The convergence is 
confirmed by ‘no change’ in the propagating mode 
amplitudes in the inflow and outflow regions. The 
waveforms are extracted from the 64 flow files 
used in the fast Fourier transform (see section 5). 
A nearly sinusoidal waveform is obtained at all 4 
points for this single frequency, 2BPF excitation. 
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BPF excitation: 

As stated above, for the fan stage (11 rotor blades 
and 27 stator vanes) considered in the current 
study, the blade passing frequency tone is cutoff. 
However, there exist two mode orders, the rotor- 
stator interaction mode m= -16 and the rotor 
locked mode m=ll, which decay exponentially 
with distance from leading and trailing edges. 
Since two modes exist, the circumferential mode 
pattern is not as clean and straight as in the case of 
2BPF excitation. Also note that the mode order 
m= -16 dominates in the inflow region while mode 
order m=ll dominates the outflow region. This 
results in the appearance of the corresponding 
mode orders in Figure 6a. The decay of these two 
modes in the inflow region (Figure 6b) and outflow 
region (Figure 6c) are shown. These cutoff modes 
decay exponentially as would be expected. 

The converged waveforms for the BPF excitation 
are shown in Figure 7. The airfoil surface pressure 
waveforms are nearly sinusoidal as for the 2BPF 
excitation. However, the inflow plane and outflow 
plane waveforms show substantial departure from 
a sinusoidal waveform. This is primarily due to 
the fact that two circumferential mode orders, 
m=l 1 and m= -16, exist for this excitation. 

Effect of Amplitude of Excitation: Single 

Frequency 

The effect of amplitude of excitation on the 
propagating mode amplitude at 2BPF is shown in 
Figures 8a and 8b, for the inflow and outflow 
regions respectively. When the excitation 
amplitude is 10 times the original amplitude 
(a 2 =0.003), the 2BPF response also increases 10 
times as shown in Figures 8a and 8b. The dashed 
line in the figures represents the response for the 
excitation amplitude for a 2 =0.003 multiplied by a 
factor of 10. The linearity of the response for the 
single frequency excitation is clearly exhibited. 

For the BPF excitation, the variation of response 
with the amplitude of excitation for the interaction 
mode m= -16 is shown in figures 8c and 8d. In the 
inflow region, the linearity is maintained for the 
range of amplitudes from 0.007 to 0.03. But in the 
outflow region, the BPF response exhibits a 
slightly nonlinear behavior which is indicated by 


the deviation of the dashed line from the solid line 
for an excitation amplitude of 0.03. 

Two-Frequency Excitation 

Next, the acoustic responses for excitation at two 
frequencies, BPF+2BPF, are considered. The 
excitation amplitudes are ai=0.03 and a 2 =0.003. 
The circumferential mode order patterns are shown 
in Figures 9a and 9b for BPF and 2BPF, 
respectively. Figure 9a for BPF, is similar to 
Figure 6a, the mode pattern observed for BPF 
alone excitation. However, the mode order pattern 
for 2BPF (Figure 9b) is distinctly different from 
that observed in Figure 4a for 2BPF alone 
excitation, in particular in the outflow region. This 
is due to the interaction between frequencies and 
self interaction as discussed below. 

The converged waveforms for the two frequency 
excitation are shown in Figure 10. Figures 10a and 
10b show the waveforms on the airfoil surface and 
10c and lOd at the inflow and outflow planes, 
respectively. Note that the waveforms are at the 
same physical locations as in the single frequency 
excitations. For two-frequency excitations, the 
waveforms show substantial departures from the 
sinusoidal form, but appear similar to that of BPF 
excitation, indicating the dominance of BFF 
excitation amplitude. 

As noted before, the BPF tone is cutoff for this fan 
stage configuration. The decaying BPF tone at 
circumferential mode orders m= -16 and m=l 1 
have the same amplitude response for single and 
two-frequency excitations in both the inflow and 
outflow regions (Figures 11 a- lid). However, the 
2BPF response is significantly different from that 
for single 2BPF frequency excitation in the inflow 
(Figure lie) and outflow (Figure Ilf). This 
marked difference appears to stem from self 
interaction and interaction between the 
frequencies. 

Effect of Amplitude of Excitation: Two-frequency 

The gust with the amplitude of BPF, ai=0.03 and 
2BPF, a 2 =0.003 is designated as Gust (see Figure 
12). The two amplitudes were decreased to 
1/2 (ai=0.015, a 2 =0.0015), 1/4 (a^O.0075, 

a 2 =0.00075), and 1/8 (a^O.00375, a 2 =0.000375), 
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successively. The responses at these excitat- 
ions are shown in Figure 12. The dashed lines 
in the figure represent the values of the response 
at the lowest (1/8) amplitude multiplied by 2, 4, 
or 8 times. For Gust amplitude 1/4 almost a 
linear response is observed. For Gust 
amplitudes 1/2 and 1, the responses deviate 
substantially from the linear multiple of the 
lowest response amplitude. This again seems to 
stem from the self interaction, as the BPF 
amplitude is increased. This is explored further 
below. 

Self interaction 

A single frequency excitation, depending on the 
amplitude of excitation, can produce higher 
harmonic responses in addition to the response 
at the frequency of excitation [4]. In the current 
study for example, high amplitude BPF 
excitation can produce response at 2BPF. This 
is examined in Figure 13, which shows the 
2BPF response for BPF excitation amplitudes 
0.015, 0.0225, and 0.03. Also shown for 
comparison is the 2BPF response for a 2BPF 
amplitude of excitation, a 2 =0.003. At the 
inflow, the 2BPF response increases with 
increase in amplitude of excitation at BPF. For 
ai=0.015 there is substantial 2BPF response 
while for ai=0.03 the 2BPF response is higher 
than that for a 2BPF excitation with a 2 =0.003! 
At the outflow, the 2BPF response increases 
with increase in BPF excitation amplitude, but 
the response amplitudes are smaller compared to 
that of 2BPF alone excitation. 

Two-Frequency Interaction 

When the excitation is at two frequencies, 
interaction between the frequencies of excitation 
can occur depending on the relative magnitudes 
of excitation [4]. The 2BPF response, for two 
frequency excitations are shown in Figure 14 for 
three cases: ai alone, a 2 alone, and ai+a 2 . Figures 
14a and 14b show the 2BPF responses when the 


amplitudes of excitations are small and equal to 
0.003, while Figures 14e and 14f show the 
response for an amplitude of 0.03. In either 
case, a nearly linear behavior is observed both at 
the inflow and outflow regions. Figures 14c and 
14d show the responses, when ai=0.03 and 
a 2 =0.003. In this case, the interactions modify 
the response amplitudes substantially. Further 
analyses are needed to quantify the interaction 
effects. 


7. CONCLUDING REMARKS 

The acoustic response of a 2-D cascade was 
examined by solving the full nonlinear Euler 
equations employing higher order accurate 
spatial differencing and time stepping methods. 
Both single frequency (BPF or 2BPF only) and 
two frequency (BPF+2BPF) excitations were 
considered to facilitate understanding of the 
results and to compare with each other. Single 
frequency excitations show a linear response 
with the amplitude of excitation for a wide range 
of excitation amplitudes. For two frequency 
excitations, significant nonlinear effects are 
observed when the BPF amplitude is high. Self 
interaction seems to contribute to this behavior. 
The results of the two frequency excitations also 
show that a linear behavior exists as long as the 
amplitudes of excitations are equal. When the 
BPF excitation amplitude is higher, interaction 
of the two frequencies and self interaction 
modify the amplitudes of the propagating 2BPF 
mode significantly, both in the inflow and 
outflow regions. To the authors’ knowledge this 
is the first time that self interaction and two- 
frequency interactions have been computed, in a 
time-domain analysis, for a realistic gust- 
cascade interaction configuration. Further 
analysis is needed to ascertain the validity of the 
observations and accuracy of the magnitude of 
the computed interactions. 
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Figure 1. — Cascade geometry definition. 
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Figure 2. — Typical cascade grid, 9 grid blocks (12,704 points) per passage, used in the simulation. 
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Figure 3. — Cascade unsteady flow features: 2BPF excitation, a 2 =0.003; 27-passage solution. 




Figure 4. — Cascade (27 passages). 2BPF gust imposed, a 2 = 0.003. Propagating mode order at 
inflow and outflow: m = -5. 
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Figure 7. — Converged waveforms for BPF excitation, ai=0.03: (a) pressure surface (mid chord), 
(b) suction surface (mid chord), (c) inflow plane (-1.5, -0.3). (d) outflow plane (1.5, 0.0). 





Figure 8. — Single frequency excitations: Effect of amplitude of excitation on propagating 2BPF 
mode and decaying modes of BPF. 
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Figure 11. — Single frequency and two frequency excitations: dashed line BPF+2BPF, ai=0.03, 
a 2 =0.003. solid line BPF or 2BPF, ai=0.03 or a 2 =0.003. 
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Figure 12. — Response amplitudes for different excitation amplitudes: Gust: ai=0.03, a 2 =0.003; 
ratios 1/2, 1/4, and 1/8 applied to both ai and a 2 excitation amplitudes. 
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